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Abstract. We analyze the Luminosity Functions (LFs) 
of a subsample of 69 clusters from the RASS-SDSS galaxy 
cluster catalog. When calculated within the cluster phys- 
ical sizes, given by r 2 oo or r 500 , all the cluster LFs appear 
to have the same shape, well fitted by a composite of two 
Schechter functions with a marked upturn and a steepen- 
ing at the faint-end. Previously reported cluster-to-cluster 
variations of the LF faint-end slope are due to the use of 
a metric cluster aperture for computing the LF of clusters 
of different masses. 

We determine the composite LF for early- and late- 
type galaxies, where the typing is based on the galaxy 
u — r colors. The late-type LF is well fitted by a single 
Schechter function with a steep slope (a = —2.0 in the r 
band, within r2oo)- The early-type LF instead cannot be 
fitted by a single Schechter function, and a composite of 
two Schechter functions is needed. The faint-end upturn 
of the global cluster LF is due to the early-type cluster 
galaxies. The shape of the bright-end tail of the early- 
type LF does not seem to depend upon the local galaxy 
density or the distance from the cluster center. The late- 
type LF shows a significant variation only very near the 
cluster center. On the other hand, the faint-end tail of the 
early-type LF shows a significant and continuous variation 
with the environment. 

We provide evidence that the process responsible for 
creating the excess population of dwarf early type galaxies 
in clusters is a threshold process that occurs when the 
density exceeds ~ 500 times the critical density of the 
Universe. 

We interpret our results in the context of the 'harass- 
ment' scenario, where faint early-type cluster galaxies are 
predicted to be the descendants of tidally-stripped late- 
type galaxies. 



1. Introduction 

The galaxy Luminosity Function (LF) is a fundamental 
tool for understanding galaxy evolution and faint galaxy 



populations. The shape of the cluster LF provides infor- 
mation on the initial formation and subsequent evolution 
of galaxies in clusters while the slope of the faint-end in- 
dicates how steeply the dwarf number counts rise as a 
function of magnitude. 

Much work has been done on the cluster LF, with var- 
ious groups finding differences in its shape and the faint- 
end slope. Different techniques have been used to measure 
LFs of individual clusters or to make a composite LF from 
individual clusters LFs (e.g. Dressier 1978; Lugger 1986, 
1989; Colless 1989; Biviano ct al. 1995; Lumsden et al. 
1997; Valotto et al. 1997; Rauzy et al. 1998; Garilli et al. 
1999; Paolillo et al. 2001; Goto et al. 2002; Yagi et al. 2002; 
Popesso et al. 2004a). Whether the LF of cluster galaxies 
is universal or not, and whether it is different from the LF 
of field galaxies are still debated issues. Several authors 
(Dressier 1978; Lumdsen et al. 1997; Valotto et al. 1997; 
Garilli et al. 1999; Goto et al. 2002; Christlein k Zabludoff 
2003) have found significant differences between the LFs 
of different clusters as well as between the LFs of cluster 
and field galaxies, while others (Lugger 1986, 1989; Col- 
less 1989; Rauzy et al. 1998; Trentham 1998; Paolillo et 
al. 2001; Andreon 2004) have concluded that the galaxy 
LF is universal in all environments. Another debated issue 
is the slope of the faint end of the LF of cluster galaxies 
(see, e.g., Driver et al. 1994; Lobo et al. 1997; Smith et 
al. 1997; Phillipps et al. 1998; Boyce et al. 2001; Beijers- 
bergen et al. 2001; Trentham et al. 2001; Sabatini et al. 
2003; Cortese et al. 2003). The LF of cluster galaxies is 
typically observed to steepen faint-ward of M g ~ —18, 
with power-law slopes a ~ —1.8 ± 0.4. This corresponds 
to the debated upturn of the cluster LF due to an excess 
of dwarf galaxies relative to the field LF. The effect may 
be real, and due to cluster environmental effects, but it 
could also be generated by systematics in the detection 
techniques of faint, low surface-brightness galaxies. 

In Popesso et al. (2004a, hereafter paper II) we have 
recently analyzed the LF of clusters from the RASS- 
SDSS (ROSAT All Sky Survey - Sloan Digital Sky Sur- 
vey) galaxy clusters survey down to —14 mag. We con- 
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eluded that the composite cluster LF is characterized by 
an upturn and a clear steepening at faint magnitudes, in 
all SDSS photometric bands. Different methods of back- 
ground subtraction were shown to lead to the same LF. 
The observed upturn of the LF at faint magnitudes was 
shown in particular not to be due to background contam- 
ination by large scale structures or multiple clusters along 
the same line of sight. We concluded that the observed 
steepening of the cluster LF is due to the presence of a 
real population of faint cluster galaxies. 

The composite LF was well fitted by the sum of two 
Schechter (1976) functions. The LF at its bright-end was 
shown to be characterized by the classical slope of —1.25 in 
all photometric bands, and a decreasing M* from the z to 
the g band. The LF at its faint-end was found to be much 
steeper than the LF at its bright-end, and characterized by 
a power-law slope —2.5 < a < —1.6. The observed upturn 
of the LF was found to occur at —16 in the g band, and 
at —18.5 in the z band. 

A steep mass function of galactic halos is a robust pre- 
diction of currently popular hierarchical clustering the- 
ories for the formation and evolution of cosmic structure 
(e.g. Kauffmann et al. 1993; Cole et al. 1994). This predic- 
tion conflicts with the flat galaxy LF measured in the field 
and in local groups, but is in agreement with the steep LF 
measured in the RASS-SDSS clusters. Two models have 
been proposed to explain the observed environmental de- 
pendence of the LF. According to Menci et al. (2002), 
merging processes are responsible for the flattening of the 
LF; the environmental dependence arises because mergers 
are more common in the field (or group) environment than 
in clusters, where they arc inhibited by the high velocity 
dispersion of galaxies. According to Tully et al. (2002), 
instead, the LF flattening is due to inhibited star forma- 
tion in dark matter halos that form late, i.e. after pho- 
toionization of the intergalactic medium has taken place. 
Since dark matter halos form earlier in higher density en- 
vironments, a dependence of the observed LF slope on the 
environment is predicted. On the other hand, if reioniza- 
tion happens very early in the Universe, this scenario may 
not work (Davies et al. 2005). Other physical processes 
are however at work in the cluster environment, such as 
ram-pressure stripping (Gunn & Gott 1972) and galaxy 
harassment (e.g. Moore et al. 1996, 1998), which are able 
to fade cluster galaxies, particularly the less massive ones. 
Whether the outcome of these processes should be a steep- 
ening or a flattening of the LF faint-end is still unclear. 

In paper II it was also shown that the bright-end of the 
LF is independent from the cluster environment, and the 
same in all clusters. On the other hand, the LF faint-end 
was found to vary from cluster to cluster. In the present 
paper (IV in the series of the RASS-SDSS galaxy cluster 
survey) we show that the previously found variations of 
the faint end of the cluster LF are due to aperture effects. 
In other words, when measured within the physical size 
of the system, given by either r2oo or rsoo, the LF is in- 



variant for all clusters, both at the bright and at the faint 
end. We also analyze how the number ratio of dwarf to 
giant galaxies in galaxy clusters depends on global cluster 
properties such as the velocity dispersion, the mass, and 
the X-ray and optical luminosities. Finally, we investigate 
the nature of the dwarf galaxies in clusters by studying 
their color distribution and suggest a possible formation 
scenario for this population. 

The paper is organized as follows. In § [21 of the paper 
we describe our dataset. In §[3]we summarize the methods 
used to calculate the individual and the composite cluster 
LFs. In §0]we summarize our methods for measuring the 
clusters characteristic radii. In §03 we analyze the result- 
ing composite and individual LFs. In §Elwe determine the 
cluster composite LF per galaxy type. In §[71 we analyse 
the environmental dependence of the LF, and compare 
the cluster and field LFs. In § [HI we provide our discus- 
sion, suggesting a possible formation scenario for the faint 
galaxy population in clusters. Finally, in §E]we draw our 
conclusions. 

For consistency with paper II and with previous works, 
we use Ho = 100 h km s _1 Mpc^ 1 , fl m = 0.3 and = 
0.7 throughout the paper. 

2. The data 

In order to study the variation of the cluster LF from sys- 
tem to system, the analysis has to be applied to a large 
statistical sample of clusters, covering the whole spectrum 
of properties (in mass, richness, X-ray luminosity and op- 
tical luminosity) of the systems considered. Since the X- 
ray observations provide a robust method of identifica- 
tion of galaxy clusters and the X-ray luminosity is a good 
estimator of the system total mass and optical luminos- 
ity (see paper I and Popesso et al. 2004c, hereafter pa- 
per III), we have used for our purpose the RASS-SDSS 
galaxy cluster sample, which is an X-ray selected sample 
of objects in a wide range of X-ray luminosity. The up- 
dated version of the RASS-SDSS galaxy cluster catalog 
comprises 130 systems detected in the RASS and in the 
SDSS sky region (16 clusters more than in the first ver- 
sion of the catalog released in paper I due to the larger 
sky area available in the SDSS DR2). The X-ray clus- 
ter properties and the cluster redshifts have been taken 
from a variety of X-ray catalogs, that allow to cover the 
whole Lx spectrum. The X-ray intermediate and bright 
clusters have been selected from three ROSAT based clus- 
ter samples: the ROSAT-ESO flux limited X-ray cluster 
sample (REFLEX, Bohringer et al. 2002), the Northern 
ROSAT All-sky cluster sample (NORAS, Bohringer et al. 
2000), the NORAS 2 cluster sample (Retzlaff 2001). The 
X-ray faint clusters and the groups have been selected 
from two catalogs of X-ray detected objects: the ASCA 
Cluster Catalog (ACC) from Horner (2001) and the Group 
Sample (GS) of Mulchaey et al. (2003). The RASS-SDSS 
galaxy cluster sample comprises only nearby systems at 
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the mean redshift of 0.1. The sample covers the entire 
range of masses and X-ray luminosities, from very low- 
mass and X-ray faint groups (10 13 M© and 10 42 ergs -1 ) 
to very massive and X-ray bright clusters (5 x 10 15 M© 
and 5 x 10 44 ergs -1 ). 

The optical photometric data are taken from the 2 nd 
data release of the SDSS (Fukugita et al. 1996, Gunn et 
al. 1998, Lupton et al. 1999, York et al. 2000, Hogg et al. 
2001, Eisenstein et al. 2001, Smith et al. 2002, Strauss et 
al. 2002, Stoughton et al. 2002, Blanton et al. 2003 and 
Abazajian et al. 2003). The SDSS consists of an imaging 
survey of ir steradians of the northern sky in the five pass- 
bands u, g, r, i, z,. The imaging data are processed with a 
photometric pipeline (PHOTO) specially written for the 
SDSS data. For each cluster we defined a photometric 
galaxy catalog as described in § 3 of Popesso et al. (2004b, 
paper I). For the analysis in this paper we only use SDSS 
Model magnitudes (see paper II for details). 

In this paper we consider a subsample of 69 clusters 
of the RASS-SDSS sample for which the masses, velocity 
dispersion, r 2 oo and rsoo (see §BJ were derived through the 
virial analysis (see paper III) applied to the spectroscopic 
galaxy members of each systems. 

Since throughout the paper the results obtained with 
the current analysis of the cluster LF are often compared 
with the results obtained in paper II, it is important to 
notice that the cluster sample used here is a subsample of 
the dataset used in paper II. 

3. Determination of the individual and composite 
Luminosity Functions 

We here summarize the methods by which we measure the 
individual and composite cluster LFs. Full details can be 
found in papers I and II. 

We consider two different approaches to the statistical 
subtraction of the galaxy background. As a first approach, 
we calculate a local background in an annulus centered on 
the X-ray cluster center with an inner radius of 3 h -1 Mpc 
and a width of 0.5 deg. 

As a second approach we derive a global background 
correction. We define as Njj g (m)dm the mean of the galaxy 
number counts determined in five different SDSS sky re- 
gions, randomly chosen, each with an area of 30 deg 2 . A 
detailed comparison of the local and global background es- 
timates can be found in paper I. The results shown in this 
paper are obtained using a global background subtraction. 

We derive the LFs of each cluster by subtracting from 
the galaxy counts measured in the cluster region, the 
field counts rescaled to the cluster area. Following previ- 
ous literature suggestions, we exclude the brightest cluster 
galaxies from the clusters LFs. 

In order to convert from apparent to absolute magni- 
tudes we use the cluster luminosity distance, correct the 
magnitudes for the Galactic extinction (obtained from the 
maps of Schlegel et al. 1998), and apply the K-correction 
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Fig. 2. The 4 panels show the composite LFs in the 4 
Sloan bands. The individual LFs used to calculate the 
composite LFs are measured within the physical sizes of 
the clusters, as given by r2oo- 



of Fukugita et al. (1995) for elliptical galaxies, which are 
likely to constitute the main cluster galaxy population. 

The composite LF is obtained following Colless (1989) 
prescriptions. A detailed description of the method can be 
found in paper II. 

3.1. Low surface brightness selection effect 

It is well known that magnitude-limited surveys may be bi- 
ased against low-surface brightness galaxies (e.g. Phillips 
& Driver 1995). An assessment of this bias for the SDSS- 
EDR and SDSS-DR1 has been done by Cross et al. (2004), 
who compared these catalogs with the Millennium Galaxy 
Catalog (Liske, Lemon, Driver et al. 2003), a deep survey 
limited in surface brightness to 26 mag arcsec -2 . Cross 
et al. (2004) concluded that the incompleteness of SDSS- 
EDR is less than 5% in the range of effective surface- 
brightness 21 < /i e < 25 mag arcsec -2 , and it is around 
10% in the range 25 < /i e < 26 mag arcsec -2 . In this 
paper, galaxies contributing to the faint-end of the clus- 
ter LFs have magnitudes 18 < r < 21. In this magnitude 
range, 65% of the objects have fi e < 23 mag arcsec -2 , 
30% have 23 < /i e < 24 mag arcsec -2 , and 5% have 
mit e > 25 mag arcsec -2 . Hence, from the results of Cross 
et al. (2004), we do not expect that the bias against low 
surface-brightness galaxies results in an incompleteness 
above ~ 5%. The faint-end of the cluster LFs derived in 
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Fig. 1. The z-band DGR vs. the cluster mass (panel a), the velocity dispersion a (panel b), the X-ray luminosity 
(panel c), and the optical luminosity (panel d). The DGR is calculated within a circle of 1 Mpc radius centered on the 
X-ray cluster center. In each panel, we list the value of the Spearman's rank correlation coefficient and the implied 
probability of no correlation. 



this paper should thus be quite unaffected by this selection 
effect. 

4. The characteristic radii of galaxy clusters 

We here describe the methods by which we measure the 
characteristic radii r 50 o and r 2 oo ■ r 200 is the- radius where 
the mass density of the system is 200 times the critical 
density of the Universe and it is considered as a robust 
measure of the virial radius of the cluster. Similarly, rsoo is 
defined setting 500 instead of 200 in the previous definition 
and it samples the central region of the cluster. Full details 
can be found in paper III. 

We estimate a cluster characteristic radius through the 
virial analysis applied on the redshifts of its member galax- 



ies. We use the redshifts provided in the SDSS spectro- 
scopic catalog to define the galaxy membership of each 
considered system. The SDSS spectroscopic sample com- 
prises all the objects observed in the Sloan r band with 
pretrosian magnitude rp < 17.77 mag and half-light sur- 
face brightness fi 50 < 24.5 magarcsec -2 . The SDSS DR2 
spectrocsopic sample used for this analysis counts more 
tha 250.000 galaxies. 

Cluster members are selected following the method of 
Girardi et al. (1993). First, among the galaxies contained 
in a circle of radius equal to the Abell radius, those with 
redshift | cz — cz c i uste r |> 4000 km s" 1 are removed, where 
Zduster is the mean cluster redshift. Then, the gapper pro- 
cedure (see also Beers et al. 1990) is used to define the clus- 
ter limits in velocity space. Galaxies outside these limits 
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Fig. 3. Contour plots of the 68%, 95%, and 99% confi- 
dence levels of the parameters of the bright-end compo- 
nent of the double-Schechter function fit to the 4 SDSS 
bands composite LFs. Solid (dotted) contours show the 
results for the composite LF computed within f2oo (re- 
spectively r 50 o). 




-18 -17.5 
Mi* -5log(h) 




-18 -17.5 
Mz* -5log(h) 



Fig. 4. Same as Fig. [31 but for the faint-end component. 



are removed. Finally, on the remaining galaxies we apply 
the interloper-removal method of Katgert et al. (2004; see 
Appendix A in that paper for more details). 

The virial analysis (see, e.g., Girardi et al. 1998) is then 
performed on the clusters with at least 10 member galax- 
ies. The velocity dispersion is computed using the biweight 
estimator (Beers et al. 1990). The virial masses are cor- 
rected for the surface-pressure term (The & White 1986), 
using a Navarro et al. (1996, 1997) mass density profile, 
with concentration parameter c = 4. This profile provides 
a good fit to the observationally determined average mass 
profile of rich clusters (see Katgert et al. 2004). 

Our clusters span a wide range in mass; since clusters 
of different masses have different concentrations (see, e.g. 
Dolag et al. 2004) we should in principle compute the clus- 
ter masses, M's, using a different concentration parame- 
ter c for each cluster. According to Dolag et al. (2004), 
c oc M -0102 . Taking c = 4 for clusters as massive as those 
analysed by Katgert et al. (2004), M ~ 2 x 10 15 M Q , Dolag 
et al.'s scaling implies our clusters span a range c ~ 3-6. 
Using c = 6 instead of c = 4 makes the mass estimates 
4% and 10% higher at, respectively, r2oo and r5oo, while 
using c = 3 makes the mass estimates lower by the same 
factors. This effect being clearly much smaller than the 
observational uncertainties, we assume the same c = 4 in 
the analysis for all clusters. 

If M V i r is the virial mass (corrected for the surface- 
pressure term) contained in a volume of radius equal to the 
clustercentric distance of the most distant cluster member 
in the sample, i.e. the aperture radius r ap , then, the radius 
r2oo is then given by: 



7*200 = r ap [p vir /{200p c )} 



1/2.4 



(1) 



where p V i r = 3M V i r / (4Trr^ p ) and p c {z) is the critical den- 
sity at redshift z in the adopted cosmology. The exponent 
in eq.Q is the one that describes the average cluster mass 
density profile near r2oo, as estimated by Katgert et al. 
(2004) for an ensemble of 59 rich clusters. Similarly, rsoo 
is estimated by setting 500 instead of 200 in eq.Q). 



5. Analysis of the individual and composite LFs 

In order to analyze the behavior of the composite LF 
faint-end as a function of waveband and clustercentric 
distance, we define the number ratio of dwarf to giant 
galaxies, DGR, as the ratio between the number of faint 
(-18 < M < -16.5) and bright (M < -20) galaxies 
in the cluster LF. The DGR is found to vary from clus- 
ter to cluster, more than expected from statistical errors. 
These variations are not random however. As shown in 
Fig.^ when the DGRs are computed within a fixed metric 
radius, they are significantly anti-correlated with several 
cluster global properties, i.e. the cluster velocity disper- 
sions, masses, and X-ray and optical luminosities (veloc- 
ity dispersions, virial masses, and X-ray luminosities for 
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our cluster sample were derived in paper III). All the cor- 
relations are very significant (1-5 x 10 -5 , according to a 
Spearman correlation test). The more massive a cluster, 
the lower its fraction of dwarf galaxies. 

The correlation between cluster DGRs and cluster 
masses is most likely due to the choice of a fixed met- 
ric aperture for all the clusters. In fact, a fixed metric 
aperture samples larger (smaller) fractions of the virialized 
regions of clusters of smaller (respectively, larger) masses, 
and DGR is known to increase with clustercentric distance 
(paper II). 

Because of this effect, the different cluster physical 
sizes must be taken into account before comparing dif- 
ferent cluster LFs. We then determine the individual and 
composite LFs within r§oQ and r2oo for the subsample of 
69 clusters of the RASS-SDSS galaxy cluster sample for 
which these parameters are known (see paper III). 

The composite LF calculated within r2oo is shown in 
Fig. [21 for four SDSS photometric bands. The it-band LF 
is not shown; in this band, there is no evidence for an 
upturn at faint magnitude levels (see paper II). For all 
the other bands LFs, a single Schechter function does not 
provide acceptable fits, and a composite of two Schechter 
functions is needed: 

m = r[{^T b exvA + (±)«texpA (2) 

where b and / label the Schechter parameters of the bright 
and faint end respectively. From the reduced-x 2 values 
given in TableQ]we conclude that a double-Schechter func- 
tion does provide adequate fits to the 4-bands composite 
LFs. Alternatively, we fit the composite LFs with a func- 
tion of this form: 

m=<t>*(^rexp(j^)ii+(±-n (3) 

In this function, (jf , L* and a are the standard Schechter 
parameters, L t is a transition luminosity between the two 
power laws and (3 is the power law slope of the very 
faint end (Loveday 1997). Both functions require the same 
number of fit parameters. However, the double Schechter 
component function provides slightly better fits than the 
Schechter+power-law function in all the Sloan bands (see 
Table UJ. 

The Double Schechter function has been used for 
the first time by Driver et al. (1994), while Thompson 
& Gregory (1993) and Biviano et al. (1995) suggested 
a Gaussian+Schechter function, to fit respectively the 
bright and the faint end of the LF. More recently, Hilkcr 
et al. (2003) used a double Schechter Function to fit the 
LF of the Fornax cluster. 

The confidence-level contours of the best-fit parame- 
ters of the bright- and faint-end Schechter components are 
shown in Figs. |3| and 0] respectively. Both results for the 
composite LF within rsoo (dotted contours) and r2oo (solid 
contours) are shown. Clearly, the best-fit Schechter func- 
tion to the LF bright-end does not change significantly 



from r5oo to r2oo (see Fig. |3J) confirming the findings of 
paper II. However, the faint-end LF steepens significantly 
(by 0.1-0.15 dex) from r 50 o to r 2 oo, and the characteristic 
magnitude correspondingly brightens by 0.3-0.4 magni- 
tudes (see Fig. EJ, thereby indicating an increasing DGR 
with radius. Our result is in agreement with the findings of 
paper II, and several other works in the literature, which 
were however mostly based on single cluster studies (e.g. 
Lobo et al. 1997; Durret et al. 2002; Mercurio et al. 2003; 
Pracy et al. 2004; see however Trentham et al. 2001, for a 
discordant result). 

While our conclusions on the composite LF agree with 
those of paper II, we find here different results concerning 
the individual cluster LFs. While in paper II we claimed 
significant LF variations from cluster to cluster, we dis- 
cover that such variations disappear when the individ- 
ual cluster LFs are computed within the physical sizes 
of each cluster (defined by r 500 or r2oo)- This can be seen 
in Fig.[SJi, where we plot the individual LFs of 15 clusters 
(those with the faintest absolute magnitude limits) and, 
superposed, the composite LF, all measured within r2oo 
and in the r-band. The agreement between the composite 
and individual LFs is very good. Fitting the composite 
LF to the individual cluster LFs result in the reduced- 
X 2 distribution shown in Fig. |3Jd. For 90% of the clusters 
the probability that the composite and individual LFs are 
drawn from the same parent distribution is larger than 
95%. 

In Fig. [SJ: we also show the z-band DGR-distribution. 
When compared to the DGR distribution found in paper 
II, the new DGR distribution is much narrower. In this 
paper we considered the DGR within r2oo of 29 clusters, 
those with known mass, r2oo and rsoo , out of the 35 sys- 
tems considered in paper II. The mean value of the DGR 
is 3.5 and its dispersion is indeed very close to the mean 
DGR statistical error of 1.4, as expected if the individual 
cluster LFs are indeed all rather similar when computed 
within a cluster-related physical radius. 

Finally, in Fig. we show DGR within r2oo as a func- 
tion of the cluster mass M200 (panel a) and the velocity 
dispersion (panel b). There is no hint of the relation pre- 
viously found (compare with Fig. [IJi) : the Spearman cor- 
relation coefficient is —0.08, which is not statistically sig- 
nificant. Similar results are found also for the DGR — Lx 
and DGR — L op relations. 

Hence we conclude that the cluster to cluster LF vari- 
ation seen in paper II arc entirely due to the use of a fixed 
metric aperture for all clusters, rather than an aperture 
sampling the same fraction of the virialized region of each 
cluster. 

6. The cluster LF per galaxy type 

In order to better understand the nature of the cluster 
galaxies responsible for the LF upturn at low luminosi- 
ties, we examine their color distribution. In particular, we 
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Fig. 5. Panel a): the individual r-band LFs within r2oo of a subsample of 15 clusters with the faintest absolute magni- 
tude limit (M r ^u m < —15.5). Empty squares and filled points distinguish the LFs computed from cluster members only 
(down to the SDSS spectroscopic completeness magnitude, r < 17.77), and using a statistical background subtraction, 
respectively. The solid line is the composite LF. Cluster names are indicated. Panel 6): the distribution of the \ 2 values 
obtained from the comparison of the composite and the 29 individual LFs of clusters with M Zi u m > —16.5 mag. Panel 
c): the z-band DGR distribution of the 29 clusters. 



use the u — r color, since the u — r distribution of Sa 
and earlier-type galaxies is well separated from the u — r 
distribution of Sb and later-type galaxies (Strateva et al. 
2001), thereby allowing to distinguish the two morpholog- 
ical samples down to very faint magnitudes. 

To define the color distribution of the cluster galaxies 
we statistically subtract the contribution of field galaxies 
(Boyce et al. 2001), using the same method applied for the 
statistical subtraction of the background from the magni- 
tude number counts. We determine the background color 
distribution of field galaxies in an annulus around the clus- 
ter with inner radius larger than r2oo ; significantly under- 
or over-dense regions (e.g. voids and background clusters) 
are excluded. By subtracting the background color distri- 
bution from the color distribution of galaxies in the cluster 
region, we obtain the u — r distribution of cluster galaxies. 
The validity of the method is confirmed by its application 
to the spectroscopic subsample, for which cluster mem- 
bership can be established from the galaxy redshifts. 



Fig.0shows the (background-subtracted) u — r distri- 
bution of cluster galaxies in the range —18 < r < —16.5 
(panel a) and —16.5 < r < —15 (panel b) for the subsam- 
ple of 15 clusters with the faintest absolute magnitude 
limit in the r band (M r ,iim > —15). The error bars shown 
in the figure take into account the galaxy counts Poisson 
statistics as well as the error due to the background sub- 
traction. 

At the redshifts of the 15 clusters considered (0.02 < 
z < 0.05) early-type galaxies have u — r colors in the 
range 2.6-2.9 (Fukugita et al. 1995), and galaxies redder 
than u — r = 3 are probably in the background. Hence, 
we can see from Fig. 0a, that the residual background con- 
tamination after the statistical background subtraction, is 
generally small (< 10 %) and in fact not significant in the 
bright magnitude range. The contamination is higher for 
the two clusters R0313 and RX 288, and probably due 
to the presence of another cluster along the same line-of- 
sight. In the fainter magnitude range, the average back- 
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Table 1. Schechter parameters of the composite LF 
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ground contamination increases to 25-35%, but is still not 
significant (see Fig. 0d). 

If we exclude galaxies with u — r > 3 from our cluster 
samples, and recalculate the cluster LFs as before (see §[3|), 
the modifications are marginal (compare filled points and 
empty squares in Fig. [SJ. If anything, a better agreement 
is now found between the composite LF and the individ- 
ual LF of the cluster R0313, for which the background 
contamination is more severe, clearly suggesting that the 
u — r > 3 color cut helps in cleaning the cluster sample 
from background contamination. 

We therefore adopt the u — r < 3 color cut to select 
cluster members, and, following Strateva et al. (2001) we 
distinguish between cluster early- and late-type galaxies 
using a color-cut u — r = 2.22. We restrict our analysis to 
the very nearby clusters (z < 0.1) to minimize the effects 
of an uncertain K-correction on the derived colors. The 
composite LFs of the early- and late-type galaxies (defined 
on the basis of their u — r colors) are shown in Fig. [5] for 
four SDSS photometric bands. The late- type galaxy LF 
is well fitted by a single Schechter function and does not 
show any evidence of an upturn at the faint end. On the 
other hand, the early-type LF looks quite different from 
the late-type LF. It shows a marked bimodal behavior 
with a pronounced upturn in the faint magnitude region. 
The best fit parameters are listed in Table [21 Such an 
upturn is then reflected in the complete (early±late) LF, 



with the late-type dwarf galaxies contributing to make the 
faint-end of the complete LF even steeper. This result is 
in agreement with Yagi et al. (2002). They determine the 
total LF of 10 clusters within 1 h Mpc radius circle. 
They find that the early-type LF cannot be fitted by a 
single Schechter function in the magnitude range from —23 
to —15, because it flattens at Mr = —18 and then rises 
again. 

7. The environmental dependence of the LFs 

In order to gain insight into the processes responsible for 
the shaping of the LF in clusters, we here examine the 
dependence of the LF on the environmental conditions. 
In particular we analyze how the LF shape, and the rel- 
ative fraction of red and blue dwarf galaxies, vary as a 
function of the clustercentric distance. Fig. 1101 shows the 
behavior of the cluster LF calculated within different clus- 
tercentric apertures, separately for the early- type (panel 
a) and late- type (panel b) galaxy populations. Distances 
are in units of T2oo- For simplicity we only plot the best fit- 
ting functions and not the data-points. The early-type LF 
is close to a Schechter function at the center of the cluster 
(within 0.2 r2oo) and shows a marked upturn afterwards. 
The location of the upturn varies from —16.2 ± 0.3 mag 
at distances < 0.3 r2oo to — 17. 4 ± 0.4 at distances < r2oo- 
The late-type LF is well fitted by a single Schechter func- 
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Fig. 6. The £-band DGR within r2oo as a function of cluster mass (panel a) and the cluster velocity dispersion (panel 
b). If DGR is calculated within r 2 oo the anti-correlation with mass (a, Lx and L op ) disappears. 
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Fig. 7. The background-subtracted u — r distribution of the galaxy members of the 15 clusters with the faintest 
absolute magnitude limit {M r ,Um > —15). a): color distribution in the magnitude range —18 < M r < —16.5; b): color 
distribution in the magnitude range —16.5 < M r < —15. 



tion at any clustercentric distance. We do not observe blue high surface-brightness objects (Blanton et al. 2004), the 

galaxies within 0.1 r2oo- Moreover, the central late- type surface brightness selection effect should be more impor- 

LF at 0.2 r 2 oo is flatter than the LFs in the outer regions tant for the late- type LF, which, once corrected, would be- 

and shows a fainter M* . Since red galaxies are mostly come steeper at the faint-end. If anything, the difference 
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Fig. 9. The composite late-type and early-type LFs in four SDSS photometric bands. The late-type (early-type) LFs 
arc displayed in the four panels on the left (respectively, right). 

Table 2. Schcchter parameters of the early and late type galaxies composite LFs 
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in slope between the faint-ends of the early- and late-type 
LFs should thus be even larger than observed. 

These results are confirmed by the analysis of the 
early- type LFs in independent clustercentric rings. We 
consider the region at distances r < 0.3 r 2 oo (the cen- 
tral ring), 0.3 < r/r2oo < 0.7 (the intermediate ring) and 
0.7 < r/r 2 oo < 1 (the outer ring). The best fitting func- 
tions of the cluster early-type LFs within these regions 
are shown in Fig. 1111 In order to emphasize the shape 
variation of the LF, all three LFs are renormalized to the 
same value. The upturn at the faint end is brighter in the 
outer ring than in the central one, confirming the previ- 
ous analysis. Moreover, the shape of the bright end of the 
cluster LF seems to be absolutely independent from the 
faint end. The values of M* and the slope of the bright 



end are consistent within the errors in the three regions 
(as found in paper II). This suggests that the process of 
formation of the bright cluster galaxies (with magnitude 
brighter that M* — 2 mag) is the same in all the cluster 
environments. Therefore, it seems unlikely that the lack 
of dwarf systems observed at the center of the cluster is 
due to a hierarchical process of formation of the bright 
central galaxies. Indeed, in that case we should observe 
also a lack of bright galaxies in the outer ring in favor of 
large amount of dwarf systems, which is not observed. 

The analysis so far provides only results about the LF 
shape. In order to quantify the relative contribution of the 
early- and late-type dwarf galaxy populations to the faint 
end of the LF, and its dependence on the environment, 
we analyse the radial (cumulative and differential) profile 
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Fig. 10. The cluster LFs within different cluster apertures in the r band per morphological type. The increment of the 
apertures is 0.1 x r 2 oo- The normalization of the fitting function is increasing at larger apertures. Panel a) shows the 
LF of the cluster red galaxy population, calculated within different clustercentric apertures expressed in unit of r2oo ■ 
Panel b) shows the same for the cluster blue galaxy population. For simplicity we only plot the best fitting functions 
and not the data points. 



of the dwarf systems in the clusters. For this, we consider 
the galaxies with — 18 < M r < —15, and to improve the 
statistics, we stack the clusters with M r ^ im > — 15 mag, 
by rescaling the clustercentric distances in units of r2oo- 
The cumulative profiles of the fractions of dwarf galaxies 
of both the early- and the late-type are shown in Fig. 112b. 
The center (< 0.4 r2oo) contains less than 30% of dwarf 
galaxies (half of them are red systems), in the selected 
magnitude range. Dwarf galaxies are more abundant in 
the cluster outskirts; the high-density environment in the 
cluster cores is hostile to dwarf galaxies. This phenomenol- 
ogy has already been observed in several individual clus- 
ters (see e.g. Lobo et al. 1997; Boyce et al. 2001; Mercurio 
et al. 2003; Dahlen et al. 2004); 

The early-type dwarf galaxies represents 35% of the 
whole dwarf population within r2oo, i-e. most of the dwarf 
galaxies are of late-type. However, the dwarf early-type 
galaxies are the dominant dwarf population region within 
0-4 r2oo, their relative fraction reaching a plateau at 
~ 0.6 r2oo, while the late-type dwarf galaxies are more 
abundant in clusters outskirts. This is confirmed also by 
the ratio between early- and late-type dwarf galaxies cal- 
culated in contiguous clustercentric rings (differential pro- 
file, see Fig. 112b). The number of early- type dwarf galax- 
ies is twice the number of late-type dwarf galaxies within 
0.2 r2oo and then decreases to 1/2 at larger distances. 



The relation between dwarf morphology and cluster- 
centric distance translates into a morphology-density re- 
lation. In Fig. U3f we show the ratio between early- and 
late-type dwarf galaxies as a function of the number den- 
sity of galaxies brighter than M r < — 18 (the bright galax- 
ies number density profile is shown in panel c of the same 
figure). As expected, the early-type dwarf galaxies dom- 
inate in high density regions, while the late-type dwarf 
galaxies are frequent in low density regions. Clearly, the 
well known morphology density relation for cluster galax- 
ies (Dressier 1980) has an extension into the dwarf regime. 

7.1. Comparison with the field 

In order to extend the morphology-density relation for 
dwarf cluster galaxies outside clusters, we extract a sub- 
sample of galaxies from the SDSS spectroscopic sample. 
We select a fairly complete sample of galaxies in the 
redshift range z < 0.02 and in the magnitude range 
-18 < M r < -16. The late-type galaxies (it - r < 2.22) 
represent the 93% of the galactic population in that range 
of magnitude, in agreement with the results of Blanton et 
al. (2004). We then calculate for each galaxy in the sam- 
ple the local density of galaxy neighbors, by counting the 
number of systems with M r < —18 mag, within 2.5 Mpc 
projected radius and ±500 km/s of the galaxy position 
and redshift. We divide the subsample in late and early- 
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Fig. 12. The fraction of red and blue dwarf galaxies as a function of the cluster environments. Panel a) shows the 
cumulative radial profile of the fraction of blue (filled points) and red (empty squared) dwarf galaxies (— 18 < M r < — 15 
mag). The fraction is defined on the total number of cluster dwarf galaxies in the considered magnitude range. Panel 
b) shows the differential radial profile of the ratio between red and blue dwarf galaxies. Panel c) shows the differential 
radial profile of the surface density of the bright cluster galaxies in clusters (M r < —18 mag). Panel d) shows the 
relation between the surface density of bright galaxies and the fraction of red and blue dwarf galaxies calculated in 
the same clustercentric ring. 



type galaxies using the color cut of Strateva et al. (2001). 
Fig. I I'M shows the number of galaxies per bin of local den- 
sity for the two galaxy types. It is clear that late-type 
galaxies (dashed dotted histogram) populate the very low 
density regions, while the early-type galaxies distribution 
(solid histogram) has a much larger spread, with 50% of 



the systems located in regions with more than 10 galaxy 
neighbors. 

It is also interesting to compare our composite cluster 
LFs with the LF of field galaxies. Blanton et al. (2004) 
have recently derived the LF of field SDSS galaxies down 
to —12 mag. Their LF have a very weak upturn, much 
shallower and at a fainter carachtcristic magnitude than 
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Fig. 8. The LFs of 4 clusters computed as for Fig. [3] 
(filled points), and by additionally excluding all galaxies 
with u — r > 3 (empty squares) . 



in our cluster LF. The faint-end slope of their LF is —1.3, 
but could be steeper (—1.5) if a correction is applied to 
account for low surface-brightness selection effects. The 
LF of blue field galaxies is even steeper, but the authors 
do not report the value of the faint-end slope. A similar 
faint-end slope (—1-5) has also been found by Madgwick et 
al. (2002) for the LF of field galaxies from the 2dF survey. 
They also noticed an upturn in the LF, due to an over- 
abundance of early-type galaxies, making it impossible to 
fit the LF adequately with a single Schechter function. A 
previous determination of the SDSS field LF was obtained 
by Nakamura et al. (2003). They found a slope of ~ —1.9 
for dlrr, consistent with the value found by Marzke et al. 
(1994) for the CfA survey. 

The faint-end slope of our late-type cluster galaxies 
LF is steeper than most field LFs for the same galaxy 
type (see Table 3 in Paper II) but consistent with those of 
Nakamura et al. (2003) and Marzke et al. (1994). Given 
the large variance of results for the field LFs, possibly 
due to the different magnitude limits adopted, or to poor 
statistics in the fainter bins of the LF (see de Lapparent 
2003 for a thorough discussion on this topic), we conclude 
there is no significant difference between the late-type LF 
in clusters and the field. 



8. Discussion 
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Fig. 11. The early- type LF calculated within three dif- 
ferent cluster regions. Only the best fitting functions are 
plotted, for simplicity, and not the data points. The LFs 
are renormalized to the same value to emphasize the shape 
variations. 



There are many observations and theoretical models in 
the literature that try to explain the formation and evo- 
lution of cluster galaxies, red dwarf galaxies in particular. 
According to the hierarchical picture for structure forma- 
tion, small dark matter haloes form before large ones. 
If one identifies the dwarf galaxies with the small dark 
matter haloes, they are predicted to origin soon after the 
structure formation began. Dwarf ellipticals would then 
be old, passively evolved galaxies. This scenario seems to 
be inconsistent with the observations of a large spread in 
age and metallicity in the clusters dwarf early-type galax- 
ies (Conselice et al. 2001,2003; Rakos et al. 2001). Hence, 
dwarf ellipticals must have had a delayed star formation 
epoch. The delay could be originated by the intense ultra- 
violet background intensity at high redshift, keeping the 
gas of the dwarf galaxies photoionized until z ~ 1, or, per- 
haps by the intra-cluster medium confinement. The intra- 
cluster medium pressure could avoid dwarf galaxies losing 
their gas content by SN ejecta. However, this possibility 
would require a much more centrally concentrated distri- 
bution of dwarf ellipticals in clusters than is observed. 

In alternative, the excess of dwarf early-type galaxies 
in clusters could origin from the evolution of field dlrr 
when they are accreted by the clusters. The evolution of 
dlrr into dwarf early-type galaxies is supported by the 
result of van Zee et al. (2004), namely that there is sig- 
nificant similarity in the scaling relations and properties 
of dlrr and dEs. A scenario where all dwarf early-type 
galaxies evolve from dlrr via disk fading does not however 
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seem possible, because many dEs in the Virgo and Fornax 
clusters are brighter than the dlrr (Conselice et al. 2001). 

Perhaps, some dwarf early-type galaxies evolve from 
dlrr and some evolve from spirals. The evolution of spirals 
into dwarf spheroidals can occur via the process of 'galaxy 
harassment' proposed by Moore et al. (1996,1998). In this 
scenario, close, rapid encounters between galaxies can lead 
to a radical transformation of a galaxy morphology. Gas 
and stars are progressively stripped out of the disk sys- 
tems, eventually leaving a spheroidal remnant, that re- 
sembles an SO galaxy or a dwarf spheroidal, depending 
on the size of the progenitor. Direct support for the ha- 
rassment scenario comes from the discoveries of disks or 
even spiral arms in dwarf early-type cluster galaxies (Jer- 
jen et al. 2000; Barazza et al. 2002; Graham et al. 2003). 
Indirect support comes from the similar velocity distribu- 
tion of dwarf cluster galaxies (Drinkwater et al. 2001) and 
gas- rich spirals and irregulars (Biviano et al. 1997), both 
suggesting infalling orbits. 

Is the harassment scenario still viable in view of our 
results? We can draw the following conclusions from our 
observational results. First, the universality of the cluster 
LF suggests that whatever shapes the cluster LF is not 
strictly dependent on the cluster properties. Second, the 
difference between the cluster and field LF seems to be 
related to an excess of dwarf early-type galaxies in clus- 
ters. Hence, there is a cluster-related process that leads 
to the formation of dwarf early-type galaxies, regardless 
of the cluster intrinsic properties. The process cannot be 
related, e.g., to the intra-cluster gas density, or the cluster 
velocity dispersion, or the cluster mass, hence, a process 
like ram-pressure would seem to be ruled out. 

The density dependence of the relative number of 
early- and late-type dwarfs suggests that the shaping of 
the cluster LF is related to the excess mean density rela- 
tive to the field, which is the same for all clusters if, as we 
have done, the cluster regions are defined within a fixed 
overdensity radius (7200 in our case). In other words, the 
transformation of spirals, and perhaps, dlrr, into dwarf 
spheroidals or dEs, seems to be a threshold process that 
occurs when the local density exceeds a given threshold. 
Judging from Fig. ^1 this threshold seems to occur at a 
clustercentric distance of ~ 0.6-0.7 r2oo- 

We have also found that the relative number of dwarf 
early- and late-type galaxies increases with decreasing 
clustercentric distance (and increasing density). Galaxies 
near the cluster center are probably an older cluster pop- 
ulation, accreted when the cluster was smaller, according 
to the hierarchical picture of cluster formation and evo- 
lution. Hence, these centrally located galaxies have had 
more time to accomplish the morphology transformation 
than galaxies located in the cluster outskirts, which are 
more recent arrivals. 

On the other hand, very near the cluster center, an ad- 
ditional process must be at work to explain our observed 
fading of the upturn of the cluster early- type LF, and the 



decrease of both the early- and the late-type dwarf-to- 
giant galaxy ratio with decreasing clustercentric distance. 
High-velocity dispersions in clusters inhibit merging pro- 
cesses (e.g. Mihos 2004), hence it is unlikely that dwarf 
galaxies merge to produce bigger galaxies at the cluster 
centers. Consistently, we find that the shape of the bright- 
end of the early-type LF does not depend on the environ- 
ment, which suggests that bright early-type galaxies are 
not a recent product of the cluster environment. In fact, 
the luminosity density profile of bright early-type galaxies 
has not evolved significantly since redshift z ~ 0.5 (Elling- 
son 2003). 

The most likely explanation for the lack of dwarf galax- 
ies near the cluster center is tidal or collisional disruption 
of the dwarf galaxies. The fate of the disrupted dwarfs 
is probably to contribute to the intra-cluster diffuse light 
(e.g. Feldmeier et al. 2004; Murante et al. 2004; Willman 
et al. 2004). 

The difference between the cluster and field LF could 
thus be explained as a difference in morphological mix, 
plus a density-dependent dwarf early- type galaxies LF, 
that, added to an invariant bright early-type LF, produces 
a more or less important and bright upturn, depending on 
the density of the environment. 

9. Conclusion 

We have presented a detailed analysis of the cluster indi- 
vidual and composite luminosity functions down to —14 
mag in all the Sloan photometric bands. All the luminos- 
ity functions are calculated within the physical size of the 
systems given by rsoo and r2oo- The main conclusions of 
our analysis are as follows: 

- We confirm that the composite LF shows a bimodal 
behavior with a marked upturn at the faint magni- 
tude range. A double Schechter component function is 
the best fit for the cluster LF. We show that calculat- 
ing the individual and the composite LF within a fixed 
aperture for all the systems introduces selection effects. 
These selection effects justify the differences observed 
in the faint end of the individual cluster LFs studied in 
paper II and the anti-correlations between DGR and 
the global cluster properties (mass, velocity dispersion, 
optical and X-ray luminosities) observed in this work. 
If the cluster LF is calculated within the physical size 
of the system (rsoo or r2oo), the differences due to aper- 
ture effects disappear and the individual cluster LF is 
well represented by the composite LF. Therefore, we 
conclude that the shape of the cluster LF is universal 
in all the magnitude ranges. 

- We use the u — r color to study the color distribution of 
the faint cluster galaxies. The color distribution con- 
firms that the contamination due to background galax- 
ies is due to field-to-field variance of the background. 
We apply the color cut at u — r = 2.22 suggested by 
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-18 < M < -16, z < 0.02 
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Fig. 13. The density distribution of neighbors of late 
(dotted histogram) and early (solid histogram) galaxies 
in the field. We select a fairly complete sample of nearby 
galaxies (z < 0.02) from the Sloan spectroscopic sample 
in the magnitude region —18 < Mr < —16. We calculate 
for each galaxy in the sample the local density of neighbor 
galaxies counting the number of systems with Mr < — 18 
mag, within 2.5 Mpc projected radius and ±500 km/s of 
the galaxy position and redshift. The sample comprises 
1561 systems. 



Strateva et al. (2001) to separate early-type from late- 
type galaxies and study the composite LF by morpho- 
logical type. We observe that the upturn at the faint 
magnitudes shown by the complete LF is due to early- 
type galaxies while the late-type LF is well represented 
by a single Schechter function. 

We study the cumulative and the differential radial 
profile of the faint early- and late-type galaxies in clus- 
ters. The faint early-type galaxies are concentrated in 
the central regions while the faint late-type galaxies 
dominate the outskirts of the systems. The analysis 
of the color-density relation in a reference sample of 
nearby galaxies selected from the SDSS spectroscopic 
sample suggests that red galaxies could be a typical 
cluster galaxy population. Our analysis show that the 
bright red population seems to have a luminosity dis- 
tribution absolutely independent from the behavior of 
the faint red galaxies in different environments. We 
observe a fading of the LF upturn toward the cluster 
core. 



- We propose to interpret our results in term of a combi- 
nation of two processes, transformation of spirals and 
dlrr into dwarf early- type galaxies via harassment, and 
disruption of dwarf galaxies near the cluster center by 
collisions and/or tidal effects. 

Whether galaxies evolve from one type to another, in 
response to the local density, to create the morphology- 
density relation, or whether the relation is established 
when the galaxies form, is still an open issue (see, e.g., 
Dressier 2004). Photometric data alone cannot provide 
conclusive indications about the nature and the origin of 
the dwarf population in cluster. In this respect, it would 
be very useful to sample the velocity distributions of a 
large set of dwarf galaxies in clusters, in order to con- 
strain their orbital characteristics as it has recently been 
done for bright cluster galaxies (Biviano & Katgert 2004) . 
If the dwarf early- type galaxies evolve from spirals, radi- 
ally elongated orbits are expected, while if dwarf early- 
type galaxies are a more pristine cluster population, their 
orbits should resemble the isotropic orbits of ellipticals. 
Additional insights may come from higher accuracy spec- 
troscopy of the dwarf galaxies, allowing to deduce infor- 
mation about their internal velocity dispersion and metal- 
licity, which could be used to put constraints on their age 
(see, e.g., Kauffmann et al. 2004; Carretero et al. 2004). 
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